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Abstract 

The homogeneous Bethe-Salpeter equation is solved in the quenched ladder 
approximation for the vector positronium states of 4-component quantum 
electrodynamics in 2 space and 1 time dimensions. Fermion propagator input 
is from a Rainbow approximation Dyson-Schwinger solution, with a broad 
range of fermion masses considered. This work is an extension of earlier 
work on the scalar spectrum of the same model. The non-relativistic limit is 
also considered via the large fermion mass limit. Classification of states via 
their transformation properties under discrete parity transformations allows 
analogies to be drawn with the meson spectrum of QCD. 
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I. INTRODUCTION 



This paper extends our previous nonperturbative studies of Quantum Electrodynamics 
in three spacetime dimensions (QED 3 ) from scalar to vector positronium states. 

Being a confining theory, the low energy behaviour of QED 3 must be dealt with non- 
perturbatively. We consider QED3 to be a simple but effective testing ground for non- 
perturbative methods commonly applied to Quantum Chromodynamics (QCD 4 ). Again our 
approach is via a solution to the homogeneous Bethe-Salpeter equation (BSE) with fermion 
propagator input from the Dyson-Schwinger equation (DSE). We consider a numerically 
solvable system of integral equations within the quenched, ladder approximation for the 
BSE and the quenched, rainbow approximation for the DSE. It is well known that the com- 
bination of ladder BSE and rainbow DSE respects Goldstone's theorem, i.e, in the chiral 
limit of zero current fermion mass, the spectrum admits a massless "pion" ||. Our four- 
component fermion version of QED 3 admits a spontaneously broken chiral-like symmetry 
leading to a doublet of Goldstone bosons. It is reasonable to expect therefore that the im- 
portant effect of chiral symmetry breaking on the light bound state spectrum will be well 
modelled by this approximation. 

We note, however, that our truncation does break local gauge covariance. Of particular 
concern is the effect of the rainbow approximation on the analytic structure of the fermion 
propagator. This approximation is known to generate ghost poles in the complex momen- 
tum plane in both QED 4 [Q and QED 3 [|]. Such poles can be an impediment to the BSE 
formulated in Euclidean space, which samples the fermion propagator over a region of the 
complex momentum plane ||. Herein we assume the position that, provided the ghost poles 
do not impinge upon this region of the complex plane, the rainbow-ladder approximation 
is a reasonable one. We acknowledge, however, that this point deserves further investiga- 
tion. Several QED 3 studies J7| have examined the effect on the spacelike behaviour of the 
fermion propagator if the fermion-photon vertex is replaced by a more sophisticated ansatz 
satisfying the Ward-Takahashi identity. These studies find the spacelike behaviour to be 
qualitatively similar to that of the rainbow approximation. Extension of these results to the 
remainder of the complex plane is beyond the scope of this current paper. We also note 
that any improvement in the fermion-photon vertex should be matched by a corresponding 
improvement in the Bethe-Salpeter (BS) kernel in such a way as to preserve Goldstone's 
theorem. That this can in principle be achieved by consistently matching a loop expansion 
of the vertex with crossed diagrams in the BS kernel has recently been demonstrated ||. 

The extension from scalar to vector states enriches considerably the spectrum of bound 
states. This facilitates an understanding of the importance of symmetry principles in deter- 
mining the dynamics of the positronium spectrum, and allows analogies to be drawn with 
the meson spectrum of QCD4. We find that there is a one to one correspondence between 
the classification of states in 4-component QED 3 in terms of "axial parity" and charge con- 
jugation on the one hand, and the conventional J PC classification of charge neutral mesons 
in QCD 4 on the other. 

We consider bare fermion mass ranging from zero to large values. For large fermion 
masses we are able to make contact with the non-relativistic limit. In this limit we extend 
our previous derivation of the Schrodinger equation as a limit of the BSE formalism to higher 
spin states, and observe spectrum degeneracies analogous to those of heavy quark effective 
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theory. 

The paper is organised as follows. In section II we look at the Bethe-Salpeter and Dyson- 
Schwinger approximations used and set out the method we employ to find the vector bound 
state masses. Transformation properties in QED 3 , with special attention given to the newly 
considered vector states, are discussed in appendix A. These transformation properties are 
necessary for understanding the structure of Bethe-Salpeter amplitudes, and classification of 
the vector bound states. The Bethe-Salpeter coupled integral equations for the vector states 
are given in appendix B. Section III describes the nonrelativistic limit for the vector states. In 
section IV numerical vector Bethe-Salpeter solutions are reported and comparisons are made 
with existing nonrelativistic limit calculations. The results are discussed and conclusions 
drawn in section V. 



II. SOLVING THE BETHE-SALPETER EQUATION FOR VECTOR STATES 

As in previous work the BS kernel is the quenched ladder approximation kernel (bare 
one-photon exchange). Again, for convenience, we use Feynman gauge and work with the 
Euclidean metric. The BSE can then be written as 

T v (p, P) = -e 2 J *fe D(p - q) lfl S( | P + q)Y y {q, P)S{- \ P + q)j„, (2.1) 

where F„(p,P) is the one fermion irreducible positronium-fermion-antifermion vertex with 
external legs amputated. The photon propagator D(p — q) in Feynman gauge is l/(p — q) 2 . 
The fermion propagator S is the solution to a truncated DSE. For bare fermion mass m this 
truncated DSE (the quenched Rainbow Approximation) is 

Z(p) = S( P r 1 -(i0 + m)=e 2 J ^Dip-q^Siq)^. (2.2) 

From here on we use the units e 2 = 1 allowable due to the fact that in the quenched 
approximation the BSE and DSE can be recast in terms of dimensionless momentum p/e 2 
and mass m/e 2 . 

We write the fermion propagator in the following general form 

S<J>)—i*vV) + '8<rt or S( P ) = t mp J + B{p2y (2.3) 
where the vector and scalar parts of the propagator are given by 

/ 2x MP 2 ) , / 2x B( P 2 ) , x 

Non-zero B signals dynamical fermion mass generation in the massless limit. Suitable ana- 
lytic fits to A and B were found in earlier work ||, namely 

A &t (p 2 ) = - 2 ai , +a 3 e-^ 2 + l, 

Bzt(p 2 ) = t^~ 2 + he- b ^ 2 + m. (2.5) 
b 2 +p 2 
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The parameters a n ,b n (functions of fermion mass) were the result of fits to iterative solutions 
to Eq. ( p. 2D and can be found in the aforementioned paper. The analytic properties of the 
propagator are summarised at the end of this section. 

In order to solve the BSE we must write it as a set of numerically tractable coupled 
integral equations. To do this, we write the bound state amplitude Y v in its most general form 
consistent with the parity and charge conjugation of the required bound state, substitute it 
into the BSE and project out the coefficient functions for the individual Dirac components. 
The general vector vertex can be found in appendix A and the integral equations in appendix 
B. 

The solution to the BSE involves iteration of the eight coupled integral equations 
Eq. ( |B.1|) . These equations can be written in the form of an eigenvalue problem. Let 

f = (a, b, c, d, e, /, g, h) T then we have 

Jdq 3 Jd\q\K(\p\,p 3 -,\q\,q 3 -M)f(\q\,q 3 -,M)=A(M) f(|p|,p 3 ;M), (2.6) 

for a given test mass M. This equation is solved for different test bound state masses until 
an eigenvalue A(M) = 1 is obtained. Such an equation exists for each symmetry case (vector 
C — — 1, vector C = +1, axi vector C = — 1 and axi vector C = +1) and for each fermion mass 
m. 

The BSE described in this section requires a fermion propagator input in the form of 
Eq. ( |2.3| ) and this needs to be available over a region in the complex momentum plane 
defined by Q 2 from Eq. (p.6j ) with q 3 and |q| real. This is the region || 



Q = \ Q 2 = X + iY 




The DSE solution should be well behaved over Q. In previous work |2| we have studied the 
analytical properties of the propagator and here we summarise briefly. 

Conjugate poles exist where the factor p 2 A 2 + B 2 appearing in the denominator of the 
fermion propagator is zero. Ref. || lists the conjugate poles arising from the fits for each 
fermion mass and the corresponding maximum bound state masses allowed. The maximum 
M allowed is the value for which the boundary of Q in Eq. (|2.7| ) coincides with the conjugate 
poles as we should not allow the singularities to enter into the BSE sampling region. If the 
singularities were to enter that region it would be necessary for compensating zeros to exist 
in the Dirac coefficients (ie: nodes in the wavefunction). Note that the region boundary ( |2.7| ) 
is bound state mass M dependent and thus one must take care as M increases during the 
bound state search. 

We have shown M that both the DSE solution and our spacelike fits Eq. (|2.5| ) have 
conjugate singularities. For small bare fermion mass m we find that the fits reproduce well 
the position of singularities. However, as m is increased beyond the scale e 2 = 1 of the 
model, the spacelike propagator functions A and B become quite level, and the propagator 
poles recede further into the timelike half of the complex plane. In this limit, numerical 
propagator fits are unable to model accurately the important analytic pole structure in the 
timelike half of the p 2 -plane. We therefore find our numerical bound state mass solutions 
contaminated by noise for m » 1. For this reason, the non-relativistic limit m — > 00 of 
the theory must be treated separately. 



4 



III. NONRELATIVISTIC LIMIT 



In previous work || a nonrelativistic limit (m — > oo) of the BSE was derived for the 
scalar states of the positronium system. In this work we follow the same procedure for the 
vector states which will allow a comparison to be made with the large fermion mass limit of 
the full BSE calculation (solution to Eq. ( B.l|) ). 

Consider first the BSE, in which we set the bound state momentum in Eq. ( |2.1| ) equal to 



= (2m + 5)iw^ where w M = (0, 0, 1) and —5 is a "binding energy". This gives (according 
to the momentum distribution in Eq. ( |2.1| )) 



(2tt 



m + - I iw + q 



T v {q) s 



m + - I iw + q 



An expansion in orders of 1/m || leads to the propagators 



S 



m + — I iw + q 



;i+7a)/2 



and 



S 



m 



iw + q 



[-(1/2)5 + iq 3 + |q| 2 /2m] + S + (g 3 , |q|) 

. (1 ~ 73)/2 

' [-(1/2)5 - tq 3 + |q| 2 /2m] + S_(g 3 , |q|) 



O 



+ o 



lnm' 



lnm 



ni 



7m- (3-1) 



(3.2) 



where, to one loop order in the self energy, 



£±(? 3 , |q|) = ~T^ ln 



|q| 

2m 



1 \ \ . 



(3.3) 



(3.4) 



In Ref. we assumed that, if the fermion self energy is calculated to all orders in rainbow 
approximation, the fermion self energy feeds back into the loop integral via the propagator 
to replace Eq. (|3.4|) by 



£±(g 3 , |q|) = -T- ln 

47T 



2m 2 H 2m 



+ s±(g 3 , |q|) 



(3.5) 



As before we assume that this equation provides an approximation to the rainbow DSE in 
the nonrelativistic limit. 

Since the vertex Y u is defined with the fermion legs truncated, and S oc |(1 ± 73), we 
find that the only relevant part of T u is the projection ^ (1 — 73) T v \ (1 + 73) and with this 
in mind, the general vector form in Eq. ( |A9| ) becomes 



I (1 - 73) 1%, P) |(1+ 73) = I (1 - 73) [ ffiGzs, |q|) u u {q) M 

+H 2 (q 3l |q|) v v (q) tf(q) 



(3.6) 



where we have introduced the linear combinations Hi = f 3 + /4/M and H 2 = f% + /7M. In 
the axivector case the general form in Eq. ([A. 10|) becomes 



(1 - 73) |(i+73)= K1-73) 



74 

75 



#i(?3, |q|) «!/(?) 

+H 2 (q 3 , |q|) ^(g)7 45 ] 



(3.7) 
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where H x = fi + ifaM and H 2 = h + if 6 M. 

Substituting Eqs. (|3.2|) , (|3.3|) and (|3.6|) into Eq. ( |3.1|) one obtains for the vector states 
two integral equations (with H l = Hi/\q\ and H 2 = —H 2 /M 2 \q\ ) 



H[ip) 



d 3 q 1 



d 3 q 1 



^(g) cos 2 fl + i^(g)sin 2 9 



-¥ 



ML 

2m 



£±(g 3> 1*3-1) 



2 • 



H[{q) sin 2 9 + H 2 (q) cos 2 



Ml 

2m 



+ E±(g 3 , |q|) 



2 ' 



(3.8) 



In these equations 9 is the angle separating the 2- vectors p and q, and H[ 2 {q) = H x 2 (?3, |q|)- 
It is clear that two independent equations can be found by applying simple linear combina- 
tions. Let F = H 1 + H 2 and G = H l 



H 2 , then 



F(p) 



d 3 q 1 
(2^F (p-qf 



F(q) 



¥ + iq 3 + l -£- + X ± (q 3 ,\q\) 



(3.9) 



This is identical to the non-relativistic limit of the axiscalar equation which has been solved 
numerically in Ref. 0, and where only a C = +1 solution could be found (where F is even 
in q.P). As we have seen, the function F in the vector case involves a linear combination 
of functions / 3 , / 4 , f 7 , and / 8 . In appendix A we see that these are all odd in q.P for the 
vector C = +1 case and even for the vector C = — 1 case. Thus it is clear that the axiscalar 
C = +1 solution is degenerate, in the nonrelativistic limit, with a vector C = — 1 state. 
We also obtain a new equation, namely 



2 (S) - 1 



G(q) 



-l6 + iq 3 + ^ + X ± (q 3 , 



|q|) 



(3.10) 



Note that the constituents of G used in the linear combination, namely / 3 , / 4 , / 7 and / 8 , 
are all odd in q.P for the C = +1 case, and even in q.P for the C — — 1 case, and so is the 
function G. 

In the axivector case we find that the surviving terms in | (1 — 73) T^ v 2 (1 + 73) can 
be collected into two linear combinations (Hi(p) and H 2 {p) defined below Eq. (|3.7j)) and 
rescaled by dividing each by |p| to find two identical decoupled equations which are also 
identical to the scalar equation of previous work. If F(p) = Hi(p)/ |p| or F(p) = H 2 (p)/ |p|, 
we have 



F(p) 



1 



p.q 



d?q 

(2tt) 3 (p-q) 2 |p| |q| 



F(q) 



¥ + iq 3 + ^ + E ± (g 3 ,|q|) 



(3.11) 



In one case F is an even function of q.P and in the other it is an odd function. The same 
equation ( |3.11| ) could only be solved with F an even function of q.P in previous work |2j 
where the scalar C — +1 solution was found. In the current case, solutions to Eq. fl3.11|) 
can be identified with the axivector C = +1 state by making the identification F cx Hi, or 
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with the axivector C = — 1 state by making the identification F oc H%. We conclude that 
the axivector C = ±1 states are mass degenerate in the m — > 00 limit. 

We should therefore see, in the large fermion mass limit, a vector state with C = —1 
degenerate with the axiscalar state with C = +1, new vector states (possibly C = — 1 and 
C — +1) which are solutions to Eq. ( |3.10| ), and axivector states with C = ±1 which are 
degenerate with the scalar C = +1 state. The degeneracies between scalar and vector states 
are analogous to those of heavy quark effective theory H], in which the hadron is insensitive 
to the heavy quark spin to leading order in the inverse quark mass. 

Eqs. Q3.9|), ( p.lL|) and ( |j.ll| ) correspond to states with orbital angular momentum £ = 0, 
2 and 1 respectively. To prove this, we show that these integral equations are equivalent to 
the Schrodinger equation of Koures in Ref. [plj. Following the working of Ref. Q], we see 
that Eqs. ( |3.9| ), ( |3.10| ) and ( |3.11D may be rewritten as in Eq. (3.27) of Ref. 0, namely 



-S 



iPl 



Hps,\p\) = 1 7^^(lp-ql) $ fe,|q|)x(p,q) (3.12) 



where x(p> q) is one of 1, or |2 (^y^j) — l|- In the final analysis the self energy XL 

serves the purpose of cancelling infrared divergence arising from the logarithmic potential 
V. 

On the other hand, consider the Schrodinger equation for a particle of orbital angular 
momentum I in (2 + 1) dimensions ||10|| : 



-5 



rn 



+ 2ReS_ \ 0(|r|)e 



,±UB 



V(r) 0(|r|)e 



±m 



(3.13) 



One easily shows that 



F.T. of 



0(|r|)e 



±m 



c/ 2 r0(|r|)e 



ir.p±ae 



<f>(\p\)e ±m > 



where we have defined ( Je is Bessel function of order 



(3.14) 



0(|p|) = 2vr/ / J,(|p||r|)0(|r|)|r|d|r| 
Jo 



(3.15) 



and 9 is the angle the vector r makes with the rj-axis, 8 P is the angle the vector p makes 
with the pi-axis. From this we find that the F.T of the r.h.s. of Eq. ( p,13|) is the convolution 
integral 



(2tt) 2 



V(|p-q|)0(|q|)^ 



(3.16) 



while the F.T of the l.h.s. of Eq. ( pTT3|) is 



5 + EL + 2ReE_U(|p|)e ± ^ 
m 



(3.17) 



Equating Eqs. (|3.17|) and (|3.16| ) we obtain 
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^y(|p-q|)0(|q|)cos^ 3 -^) 

^nip-q|)0(|qlMp,q). (3.18) 



which is equivalent to Eq. (|3.12|) . 



To summarise, we have obtained equations for degenerate t = vector and axiscalar 
states (Eq. ( |3.9|) ), a new £ = 2 vector state (Eq. (|3.10|) ), and degenerate axivector and scalar 
states with i = 1 (Eq. ( p. 11 )). In order to explain the occurrence of this particular set of 



angular momentum and parity combinations we need to study the symmetry properties of 
bound states in (2 + 1) dimensions. 

To classify the bound states in terms of the spin and orbital angular momentum parts, 
we consider briefly the finite dimensional representations of the (2 + 1) dimensional Lorentz 



group SO{2 } 1) [11]. Its Lie algebra is defined by the commutation relations 

=ie M /£ p , (3.19) 

where S is the generator of rotations in the X\-Xi plane and S 12 generate translations in the 
Xi ; 2 directions. The finite dimensional representations are classified by the Casimir operator 

£ 2 = (£o) 2 -(£i) 2 -(£ 2 ) 2 , (3.20) 

which takes eigenvalues L(L + 1), while S takes eigenvalues i = —L, —L + 1, . . . , L — 1, L. 
Identifying £ with the angular dependence of the Schrodinger wavefunction discussed above, 
we shall refer to i as the total orbital angular momentum. To the spin part of the bound 
state angular momentum we also attach a quantum number S, which we call the total spin. 
For the symmetric spin combinations of two Dirac particles we have S = 1, and for the 
antisymmetric combination, S = 0. The spin angular momentum in the x-y plane takes 
values m = —S, —S + 1, . . . , S — 1, S. 

The bound states are also classified by their eigenstates under the action of the "axial 
parity" operator A: 

A: x 11 = (x°,x\x 2 ) x'» = (x°,-x\-x 2 ), (3.21) 

which is distinct from the parity operator defining the operation Eq. ( |A.1|) in appendix A. 
Eq. ( |3.21| ) is a rotation through n, and is the operator relevant to interchange of constituent 
fermions. The transformation properties of Dirac spinors and bilinear currents under axial 
parity are given by Eqs. ( |A.(j| ) and ( |A.7] ). The axial parity A of each of the bound states 
considered herein is listed in Table [IJ 

In the non-relativistic limit, the generalised Pauli exclusion principle [F2| adapts from 



(3+1) to (2+1) dimensions to give the following relations between the total orbital angular 
momentum L, the total spin S, axial parity A and charge parity C : 



A={-l) L+ \ (3.22) 
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C = 



(3.23) 



Two fermions can form a system with a total spin of S = or S = 1. For a specific total 
angular momentum J, vector addition of angular momenta m + I determines the possible 
orbital angular momenta £. Table | lists all possible scalar and vector states. 

Consider the scalar states where J = 0. For an orbital angular momentum of zero (L = 0) 
the total spin must also be zero (S = 0). In this case the axial parity must be negative 
(A = —1) and charge parity positive (C = +1). We know from our integral equations that 
this state is the axiscalar state. It is also known from the transformation properties outlined 
in appendix A that the axiscalar does indeed have A = — 1. The other possibility with J = 
is the L = 1, S = 1 state where the axial parity is positive and the charge parity is positive. 
This is the scalar state which we know from its transformation properties has A = +1. So 
we have seen that it is only possible for scalar states to have L = or L = 1 and they must 
have positive charge parity. Negative charge parity scalars are forbidden by the generalised 
Pauli Exclusion principle. States such as the J AC = state, for example, we refer to as 
having unnatural parity. 

Now consider the vector case (J = 1). If L = then S must equal 1. In this case both 
the axial parity and charge parity are negative. We know that this is the L = vector state. 
With no spin-orbit coupling contribution in our e~-e + potential this state is degenerate with 
the L = axiscalar state. 

If L = 1 then there are two possibilities. A J = 1 is possible with either S = or S = 1. 
Both of these cases correspond to the L = 1 axivector states, both with positive axial parity, 
one with negative charge parity and the other with positive charge parity. Again, with no 
spin-orbit coupling the two are degenerate. They are also degenerate with the L = 1 scalar 
state. 

The only other possible state is the L = 2, S = 1 state which has negative axial parity 
and negative charge parity. This is the L = 2 vector state found earlier in this section. 
There appears to be no possibility of an L = 2 vector state with positive charge parity and 
so we would hope that there is no such solution in our nonrelativistic BSE calculations. 

In general, we find that the axial parities determined from the angular momenta and 
Eq. (|3.22| ) match up with those found by analysing the transformation properties of bilinear 
currents in Appendix A. We also note that the "axi" transformation in (2 + 1) dimensions 
parallels the conventional parity transformation in (3 + 1) dimensions, and that a one to 
one correspondence between the allowed J AC states in (2 + 1) dimensions and J PC states in 
(3 + 1) dimensions follows from the generalised Pauli exclusion principle. 



IV. NUMERICAL RESULTS 

In this section we report our numerical solutions to the Bethe-Salpeter coupled integral 
equations of Eq. (B.l| ). Comparisons with the nonrelativistic calculations of Refs. |^3| , |10| 
will be made for large fermion mass. 

As in the scalar calculations of Ref. 0, a grid of 25 x 25 (|q|,<?3) tiles was used for 
the iterative solution to the vector equations Eq. ( |B.1| ) with the use of linear interpolation 
on each of those tiles for the sums (r rl , r r2 ) which are supplied at the corners of the tiles 
from the previous iteration. Again the tiles were non-uniform in size and an upper limit 
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to the momentum components (|q| and q^) were made large enough so that results were 
independent of their values. For each symmetry and a range of fermion masses (0 — 5) the 
equations were iterated leaving a set of bound state masses. 

Table [H] shows the bound state masses for each of the four nondegenerate states which 
are the vector C — +1, vector C — — 1, axivector C = +1 and the axivector C — — 1 states. 
We employ the notation J AC introduced in the last section to classify these states. The 
degeneracy under space reflection parity in 3-dimensions is assumed from this point on and, 
for example, the vector state refers to the doublet vector / pseudovector. Fig. 1 displays the 
solutions M for fermion mass 0-0.1. Fig. 2 shows M — 2m over the greater range of 0-1. 
We note that the poles in the fermion propagator fits lie outside the BS integration region 
Q for all solutions obtained. This has been the case in both scalar and vector calculations. 

For small m the bound state masses rise rapidly with increasing fermion mass. The 
J AC = 1 state is the lowest energy state followed by the 1~ + and then the 1 ++ and V 
states respectively. At some point between m = 0.049 and 0.064 the 1 _+ and the 1 ++ states 
cross. The curves then level off for masses approaching the nonrelativistic limit m > 1. 

We found that the two axivector states were near degenerate for small m. However, 
numerical problems prevented negative charge parity solutions to the axivector equations 
for m = 0.036 and higher. The eigenvalue A(M) in Eq. ( |2.6| ) splits into complex conjugate 
pairs past this value of m. We see no theoretical reason for the disappearance of the negative 
charge parity axivector state and we attribute it to deficiencies of the bare vertex, ladder 
approximation. 

As expected, for larger fermion masses the bound state mass rises predominantly as twice 
the fermion mass plus possible logarithmic corrections. However, there appears to be a good 
deal of noise in the large m solutions, reflecting the difficulty in accurately modelling the 
fermion propagator deep into the timelike region from spacelike fits. It is for this reason 
that Figure 2 does not display the larger fermion masses. Despite this noise, the ordering of 
the states remains unchanged as m becomes large. 

It is clear that this approach to the large fermion mass end of the spectrum is not suitable 
and that a special nonrelativistic treatment is required. What is required is solution to 
Eqs. ( |3.9| ), (|3.10| ) and ( |3.11| ). However, this involves the solution for the rainbow self 



energy S + from a DSE formulated in the large fermion mass limit. This is not within 
the scope of this paper and is the subject of future work. In previous work for the scalar 
positronium solutions a 1-loop approximation to this self energy was employed. This 
appeared warranted because the approximation gave a good fit to the rainbow propagator 
for spacelike momenta. However, it was found that the 1-loop approximation resulted in 
bound state masses that were out by a logarithmic correction (because of the noncancelling 
infrared divergences) and as a result the masses could not be compared to other large 
fermion mass solutions. In this work, we do not attempt to match 1-loop nonrelativistic 
solutions with our large m relativistic solutions. Despite this, we believe that the existence 
or nonexistence of a solution to the 1-loop nonrelativistic BSE equations is evidence of the 
existence or nonexistence of those positronium bound states. Therefore, Eqs. ( |3.9| ), ( |3.10| ) 
and (|3.11|) were solved with S + replaced by the 1-loop result in Eq. (|3]4|) neglecting terms 



of order -K. 

Because of the degeneracies pointed out in section III, solutions exist to Eqs. ( |3.9|) and 
( |3. 11| ) and from this we could conclude that the vector £ = 0(1 ) and the axivector £ = 1 
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^1 ± ) bound states do exist. We were left to solve Eq. ( |3.10|) for possible positive and 



negative charge parity £ = 2 vector solutions. It was found that the 1 solution did exist 
and had a mass increasing linearly in lnm as did the smaller £ solutions. However, the 1~ + 
solution mass did not level off to the same slope and appeared to have a strongly divergent 
mass in the nonrelativistic limit. We interpret this by saying that there is no 1~ + bound 
state, in agreement with the generalised Pauli exclusion principle. 

Although the 1-loop nonrelativistic masses cannot be quantitatively compared with 
masses found with the rainbow propagator, we again point out that the results do contain 
useful qualitative information. We find that the ordering of the states in the nonrelativistic 
limit is 1 (£ = 0) followed by l +± (i = 1) and then 1 (£ = 2). Based on both the rela- 
tivists calculations and the nonrelativistic 1-loop exercise, it appears that it is energetically 
favourable for the scalars to have positive charge parity, and for the vectors to have negative 
charge parity, in agreement with the QCD 4 meson analogy. 

We compare our relativistic integral equation results with those of existing QED 3 
Schrodinger equation results. Tarn, Hamer and Yung perform an analysis of QED 3 



from the point of view of discrete light cone quantisation. In the non- relativistic limit, their 
bound state masses M are solution to the Schrodinger equation 

(-—V 2 + — (C + ln(mr))l 0(r) = (M — 2mWr), (4.1) 
I m 2ir ) 

where C is Euler's constant and m is the bare fermion mass. This differential equation is 
solved for the completely symmetric states and the masses are given by the expression 

M = 2m + — lnm+ — (x' ln-V (4.2) 

4tt 2tt V 2 ir) K J 



Koures |nj solves the same equation but also for nonzero angular momentum £. For < 



< 4 this reference provides the lowest five eigenvalues (A) which, after the transformation 

A' = A + ln(2) + C (4.3) 



can be used in Eq. ( |4.2|) above. The eigenvalues to be used here are A' = 1.7969 (£ — 0), 
2.6566 {£ = 1), 2.9316 (£ = 0) and 3.1148 (£ = 2). 

The lowest £ = 0, 1 and 2 and the first excited £ = Schrodinger equation results of 
Refs. [[yj, |1(J are also plotted in Fig. 2. One could associate the £ = curve with the 



vector 1 state. The £ = 1 curve can be seen to match up with one or more of 1~ + , 1 ++ 
or l + ~. We know from our nonrelativistic analysis that the l +± (axivector) states are the 
only £ = 1 states and become degenerate in the large m limit. It therefore seems as though 
there is a surprisingly good match up between relativistic and nonrelativistic £ = 1 states. 
The lowest £ = 2 and the first excited £ = curves should coincide with higher vector 1 
states but no such states were solved for in our relativistic BSE treatment. Also there is no 
nonrelativistic 1 _+ state to match up with the relativistic vector state with positive charge 
parity as shown in our nonrelativistic BSE exercise and as predicted by the generalised 
Pauli exclusion principle. We emphasise that the Schrodinger equation results are of course 
formulated in the nonrelativistic limit and we cannot expect close numerical agreement with 
the small mass relativistic results. In fact, there is quite a broad region of intermediate 
masses where neither the relativistic BSE nor the Schrodinger results are suitable. 
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Finally, we may make qualitative comparisons between the calculated spectrum of J AC 
states in QED3 and the spectrum of observed J PC meson states in QCD4. Since the observed 
meson spectrum is best established for light mesons, the comparison is made in the chiral 
limit m = 0. In Table |T| we list the m = spectrum of scalar and axiscalar states from 
Ref . [0] , and of vector and axivector states from this work, together with the corresponding 



observed light meson states f!4 |. (Note that, as pointed out in the previous section, the 
axial parity quantum number A in (2 + 1) dimensions plays the role of conventional parity 
P in (3 + 1) dimensions.) Firstly we see that the gaps in the observed meson spectrum 
occur for parity combinations disallowed by the generalised Pauli exclusion principle. In our 
fully relativistic QED 3 calculations, there is nothing in principle to prevent such states from 
occurring away from the nonrelativistic limit. We are unable to say whether the unnatural 
parity solutions are a genuine property of QED 3 , or an artefact of the rainbow, ladder 
approximation. If these states are ignored, we observe surprising agreement between the 
ordering and relative magnitudes of bound state masses in QED 3 and the observed meson 
spectrum. In both cases the dynamics of the bound state spectrum is driven mainly by 
confinement and chiral symmetry breaking, these being the common features of QED3 and 
QCD 4 . 



V. CONCLUSIONS 

We have extended our earlier study of the positronium states of 4-component QED3 
from scalar to vector states. QED3 is a confining theory, and the positronium states are 
in some sense the analogues of mesons in QCD 4 . In (2 + 1) dimensions the bound states 
are classified as eigenstates of reflection parity P, axial parity A (i.e. a spatial rotation 
through 180°), and charge conjugation C. In the four component version of massless QED3, 
a £7(2) symmetry analogous to chiral symmetry is spontaneously broken to £7(1) X £7(1). The 
resulting spectrum consists of reflection parity doublets for which P is an exact symmetry. 
The states are therefore classified by the quantum numbers J AC '. We had previously studied 
in detail axiscalar (0~) and scalar (0 + ) positronium. Herein our focus is mainly on the 
extension of this work to vector (1~) and axivector (1 + ) states. 

We calculate bound state masses using the combination of rainbow Dyson-Schwinger and 
ladder Bethe-Salpeter equations. Fermion propagators calculated for spacelike momenta 
from the Dyson-Schwinger equation are extended into the required part of the complex 
momentum plane by making analytic fits to the spacelike part of the propagator. There are 
two important issues raised by this extrapolation procedure, both of which we addressed for 
the scalar case in our previous work. 

Firstly, the solution to the rainbow approximation Dyson-Schwinger equation, and the 
analytic fits, contain ghost poles in the timelike half of the complex momentum plane. It 
is necessary for the success of the approximation that these poles do not impinge on the 
set of complex momentum values sampled by the BSE. As for the scalar case considered 
previously, we find that the vector states do not pose a problem in this regard. This is 
because the scalar and vector bound state masses are typically of comparable size over the 
range of bare fermion masses considered. 

Secondly, in both the scalar and vector cases, the extrapolation procedure is inadequate 
for intermediate or large bare fermion masses (viz. m/e 2 > 0.5). This is because the 
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important contributions to the BSE in the heavy fermion limit come from deep into the 
timelike part of the complex plane. For this reason the heavy fermion limit must be treated 
separately. 

We have obtained numerical solutions to the combination of Dyson- Schwinger and Bethe- 
Salpeter equations over a range of dimensionless bare fermion masses m/e 2 for each of the 
axial and charge parity combinations 1 , 1~ + , l + ~ and 1 ++ . Solutions exist over the broad 
range of bare fermion masses considered except for the l +_ , for which the eigenvalue of the 
integral operator in the BSE becomes complex for m/e 2 > 0.036. We interpret this as a 
shortcoming of the rainbow-ladder approximation. 

In our previous work we obtained the nonrelativistic (i.e. m —>■ oo) limit of the BSE for 
scalar and axiscalar states in the form of a Schrodinger equation. Here we have extended the 
proof to vector and axivector states in a way which enables us to identify the orbital angular 
momentum of each state. We find that the relation A = (— 1) L+1 is automatically satisfied, 
and are furthermore able to restrict the allowed charge parities by using the generalised 
Pauli exclusion principle. In this limit it becomes clear that it is axial parity, and not space 
reflection parity, which is the QED3 counterpart of conventional parity in QCD4. 

We also find that in the nonrelativistic limit, the bound state spectrum becomes insen- 
sitive to the spin of constituent fermions. This is manifested as a degeneracy between 0~ + 
and 1 states and between the ++ and l +± states. Precisely the same phenomenon occurs 
in the observed heavy meson spectrum to within order of the inverse heavy quark mass, and 
can be explained in terms of heavy quark effective theory || . 

While certain J AC combinations are disallowed by the generalised Pauli exclusion prin- 
ciple in the nonrelativistic limit, there is nothing in principle to prevent their occurrence 
in the relativistic regime. We have in fact obtained unnatural parity solutions in both our 
earlier scalar and current vector positronium calculations. This is not in agreement with the 
observed meson spectrum, even for light quark mesons, where for instance negative charge 
parity scalar or pseudoscalar mesons are never observed. One is led to question whether the 
unnatural parity solutions in QED3 are an artefact of the rainbow-ladder approximation. 
More importantly we find that, if the unnatural parity solutions are ignored, there is sur- 
prising agreement between the relative mass scales of calculated J AC states in QED 3 and 
observed J PC mesons in QCD 4 , as evidenced in Table 

The original aim of studying 4-component QED3 was to explore nonperturbatively a the- 
ory which has properties in common with QCD4, namely confinement and chiral symmetry 
breaking, but without the complications of being nonabelian. We have demonstrated that 
the Bethe-Salpeter formalism generates a bound state spectrum with qualitative characteris- 
tics in common with the observed meson spectrum. One shortcoming of the rainbow-ladder 
formalism is that gauge covariance is broken by using a bare fermion-photon vertex. In any 
practical truncation it is likely that some symmetry of the original problem will be lost. In 
this paper we have adopted the position that it is more important for bound state mass 
calculations to maintain the chiral symmetry breaking mechanism than to maintain gauge 
symmetry. An alternative approach is that of lattice gauge theory, in which gauge symmetry 
is maintained, but the original chiral symmetry of QED3 survives only as a small remnant 
subgroup J0|. Specifically, lattice gauge theory simulations of quenched, non-compact QED3 
with 4-component fermions can be carried out with comparative ease fL6 |, and the chiral 



condensate checked against the predictions of model Dyson- Schwinger calculations [O. We 
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are unaware of any existing bound state spectrum calculation of this simple lattice model. 
Such a calculation would provide a much needed cross check between two vastly differing 
non-perturbative field theory methods. 



APPENDIX A - GENERAL VECTOR AMPLITUDE IN QED 3 

Solution to the BSE in Eq. fl2.1| ) for the vector states requires a general vertex function 
which will allow individual Dirac components to be projected out leaving a system of coupled 
integral equations. We assume a knowledge of the transformation properties of QED 3 
and of the Scalar BS amplitudes ||. 

We work in the four-component version of QED3 where there exists a complete set 
of 16 matrices (where [i = 0, 1 and 2) {7^} = 74, 75, 745, 7« 7u4> 7*5) 7*i5} satisfying 
^tr(7 J 47' B ) = 5 a- These matrices can be found in Ref. [0J. The QED 3 action is invariant 
with respect to discrete parity and charge conjugation symmetries, which for the fermion 
fields are given by 

i){x) -> ifj'(x') = Uip(x), ii{x) -> ^V) = ^(x)rr\ (A.i) 



il){x) -»■ ip'(x) = C?p(x) T , ip{x) -> ip'(x) = -ip(x) T C\ (A.2) 

where 2 ). The parity and charge conjugation matrices (II and C respectively) 

are each determined up to an arbitrary phase by the condition that the action be invariant [|TJ : 

II = 7l4 e^ 745 , C = 7 2 e i0c745 , (0 < <j> P , 4> c < 2tt) ( A.3) 

Vector, pseudovector, axivector and axipseudovector bound states are defined by the 
following transformation properties under parity transformations 

K pv ^) - K pv '^') = -Rpa;^ pv (x), (a.4) 

Where 

R P =(- COs2 9 t P ~ S[n ^ P ) (A-5) 

sm20p cos20p J y ' 

and in Minkowski space = diag(l, — 1, 1). Similar transformation properties exist for 
charge conjugation. 

We also find the need to classify our nonrelativistic states in terms of the transformation 
properties under what we call "axial parity" which is the (2+l)d analogue of the (3+l)d 
parity where ). The fermion fields transform like 

ij){x) -> W) = S^(x), $(x) -> il{x') = ^(x)S~ l , (A.6) 
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where S n is the operator which corresponds to a rotation through angle 7r in the x 1 x 2 plane. 
We find that a suitable operator S n which performs such a rotation and leaves the action 
invariant is the matrix Z7 . This is the same operator as used in the (3+l)d case. The 
phase i is responsible for making S% = — 1 and has no effect on the transformation of 
bilinear currents of concern here. Given that the bound state wavefunctions transform like 
the bilinear currents J a = ipixjlAipix) we find the following transformation properties for 
scalars and vectors under axial parity. 



$ s (x) - 


-> ® s (x), 






$ P5 (x) - 


-> $ P5 (x), 






$ A5 (x) - 


-> - (^(x), 


$f (x) - 


- -T%*i V (x), 


§ APS {x) - 




$f^(x)- 





(A.7) 

The matrix rj^ is defined as diag(l, —1, —1). The corresponding axial parities (A = ±1) 
for each state are listed in Table |. This definition parallels the (3+l)d case where a vector 
state transforms like rf $ v and is assigned negative parity. 

Let us now derive a general vector vertex. Consider the momentum p^ which transforms 



like p M — > A^p„ under parity. The bound state momentum P M seen in Eq. (2A) is another 
such vector. Three other vectors can be found which transform in this way which are 7^, 
u^ip) and v At (p)745- The vectors and v{p) used in previous work [|I|] are mutually 
orthogonal with P and are defined as 

P 2 p-(Pp)P Pxp 
U{P) = PV - (RpJ ' (A ' 8) 

The proof that these vectors transform like p is simple using the knowledge that a vector 
V^ip, P) transforms under parity like V^(p, P) — > UV^(Ap, AP)IT _1 and that the five matrices 
7^, 74 and 75 form an anticommuting set and 745 = —27475. 

So we have a set of five vectors with even parity. A linear combination of P and u can be 
found which gives p and so it is not required. The four vectors we work with are thus P M , 7^, 
u^{p) and v M (p)745. Any other available vectors even in parity would be linear combinations 
of these. 

A massive boson of spin 1 will have a purely transverse wavefunction Thus the 

vector vertex will be purely transverse. This zero divergence criterion in momentum space 
means that P^Y M = 0. We must find linear combinations of the four vectors above which 
are all purely transverse. Only three can be found which are T M (P) = 7 M — P M P/P 2 , u^(p) 
and v M (p)745- It therefore appears as though the most general vector T^p, P) we can write 
which is even under parity and purely transverse is a sum of these quantities, each multiplied 
by the most general scalar vertex ||J. The result is an expression with twelve terms and 
with twelve parameters. However, some terms have been included more than once. It can 
be shown that four of the twelve terms can be eliminated leaving the general vector vertex 

rl(p, p) = «» (/1 + h r + fs m + h Mi*) 

+ « m (p)745 (is + UP + h Hp) + h Ap)i^) • (A.9) 

where f n (n =1-8) are functions of q 2 , P 2 and q ■ P only. The pseudovector (PV), axivector 
(AV) and axipseudovector (APV) vertices are given by 



15 



r7(p^) = 74 5 r^( P ,p), 

rf y (p,P)= 7 45 r^(p,P) = z(~] 5 )r>,p). (A.io) 

For a specified charge parity C — ±1 of a bound state, the parity of the Dirac coefficients 
(/„) under the transformation q ■ P — > — g • P can be determined. The quantity q ■ P is the 
only Lorentz invariant which changes sign under charge conjugation and thus determines 
the charge parity of those functions. Table [IV] lists the charge parities of each function for 
each case. 

We can see from this table that multiplying the vector vertex by 745 (and thus producing 
the pseudovector vertex) leaves the function charge parities unchanged but multiplying the 
axivector vertex by 745 (producing the axipseudovector vertex) reverses the function charge 
parities. We know that the same BSE results when we multiply the vertex by 745. This 
means that the vector C — +1 and pseudovector C — +1 states are degenerate as are the 
vector C — — 1 and pseudovector C — — 1 states. We must also find the degenerate pairs 
axivector C — +1 / axipseudovector C — —I and axivector C = —1 / axipseudovector 
C = +l. 

Our conventions for Euclidean space quantities are summarised in Appendix A of Ref . [Q . 
In particular Euclidean momenta and Dirac matrices are defined by 

p(E) _ . p (M) p (E) _ p (M) (E) _ (M) (E) _ . (M) 

r 3 ~ "0 > -'1,2 _ r l,2 1 73 — 70 J 7l,2 _ '71,2 ■ 



APPENDIX B - COUPLED INTEGRAL EQUATIONS 



Working in the rest frame of the bound state where P M = (0,0,iM), the vector and 
axivector vertices given in appendix A are used to derive coupled BS equations. After 
considerable work the eight coupled equations were found to be (after rescaling and angular 
integration) 



a(p) 

Hp) 

c(p) 

dip) 
e(p) 

f(p) 

g(p) 



(2tt) 2 
1 



dq 3 



-oo 
oo 



(27r) 2 J-oo 
1 



(27r) 2 J-oo 
1 



(27T) 2 

3 

i 

(2tt) 2 
1 



dq 3 



dq 3 



dq 3 



-oo 

DO 



dq 3 



-oo 
oo 



dq 3 



(2vr) 2 J-oo 



-oo 
oo 



dq 3 



|q 
|q 
|q 
|q 
|q 
|q 
|q 



X(a,P) r al {q) 
X(a,p) r bl (q) 

[Z(a,P) r cl (q)+Y(a,(3) r c2 (q)} 

[Z{a,f3)r dl (q) + Y{a,f3)r d2 (q)) 
X(a,(3) r e2 (q) 
X(a,f3) r f2 (q) 

[Y(a,P)r gl (q) + Z(a,P)r g2 (q)} 
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h{p) = ^fi-J q " Jo |q|d|q| [ Y ( a ^) T hx(q) + Z(a,^r K2 (q)] (B.l) 

where we have rescaled the functions fx, f 2 , / 3 , ik, fa, f&, f-j and / 8 seen in appendix A to the 
new functions a, b, c, d, e, /, g and h to ensure all quantities are real. The rescalings for the 
vector case were 

/5O3, IpI), 



o(P3, IpI) = 


/i(p3, |p|), 


e (p3, |p|) 




Aba, IpI), 


/(ps, IpI) 


c(p 3 , IpI) = Ai 


'/s(p3, IpI), 


0(ps, IpI) 


%3,| P |) = M|p| 


7/4(^3, IpI), 


%3, IpI) 



M|p 

1 „ 



M| P | 2 



Mps, IpI), 



M 2 |p 



/ 8 (p 3> |p|), (B.2) 



and for the axivector case 



2 

a(p3,|p|) = jprAOs, IpI), e(p 3 , |p|) = j^Mps, IpI), 
&(ps, IpI) = ^ / 2 (p 3 , IpI), /(ps, IpI) = p /•(**> IpI). 
c(?3, IpI) = 7^/3^3, IpI), g(pa, |p|) = — ; , 2 iVfe, IpI), 

IPI M|p| 
d (PzAp\) = U(P3,\P\), %3, |p|) = 72 fs(P3, |p|)- (B.3) 

Af |P| M^Ipi 

These functions also occur within the sums (r) seen in Eq. ( p.l| ) which are defined as (for 
a function r = a-K) 

T rl (q) = T ra a(q) + T rb b(q) + T rc c(q) + T rd d(g) and 

r r2 (g) = T re e(g) + T rf f(q) + T rg g(q) + T rh h(q). (B.4) 



Also 



A>, (3) = 1 - „ P ; , , y (a, 0) 



Z(a,/3) = 1 , (B.5) 

a(a 2 -/3 2 )^ + (a 2 -/3 2 ) 



where 



a = (pa - 9s) 2 + IpI + |q| , /9 = -2|p||q|, 
The momentum Q is defined by 



Q 2 = ?3 + |q| 2 -^ 2 + ^3, (b.6) 
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and we use the abbreviations oy = oy(Q 2 ) and as = cr,s{Q 2 ) for use in the definition of the 



functions T aa , T ab , . . . which are a 


nalytic functions of q 3 , q , and M. lhe diagonal 1 s are 


given by 








T aa = (-M 2 


+ gf + |q| 2 )M 2 


TWs\ 2 , T ee = (iM 2 + g| + |q| 2 )| C r y | 2 ±| ( T 5 | 2 


T bb = — (^M 2 


+ 13 - 


-|q| 2 )Kf 


± \a s \ 2 , T ff = -{\M 2 + qi - \q\ 2 )\a v \ 2 ± |<j s | 2 


T cc = {\m 2 


+ 13- 


-|q| 2 )Kf 


±|<7 5 | 2 , T gg = (lM 2 + g 2 -|q| 2 )| ( T y | 2 ±| ( 7 5 | 2 


T dd = (\m 2 


+ gf + lql 2 )M 2 


±\a s \ 2 , T hh = {\M 2 + q 2 + \^)\a v \ 2 ±\a s \ 2 (B.7) 


where the upper sign applies to the vector equations and the lower sign to the axivector 


equations. The nonzero off diagonal T's for the vector case are 


Tab 




— T ba 


% 

= 2^ a v a s - cr* s a v )M - (a v a s + <J* s a v )q 3 


T 

J- ac 




- T 

ca 


— \H\\ a V a S + a S a V ) 


T a d 




= —Tda 


= M\q\\a v \ 2 


Tbc 




= T cb 


= 2q 3 h\Wv\ 2 


Tbd 




= T d b 


= i\q\(a v a s - cr*s a v) 


Ted 




= Tdc 


— 7:( a v a s + a%av)M + iiatvs — 0'%( T v)l3 


T 






]\/T \ rt \ |_ |2 

— 1V1 |q| |(7y 1 








— i\ ( l\( IJ v a s — °~s°~v) 


T 

J- eg 




= T 


1 / jfc \7l/T * / ^ ^ \ 

= -^v^s + <7 s <7y)M - «(eryo- s - a s a v )q 3 


T c h 


— T hc 


ol |2| 1 


1 

2 =-(^M 2 + g 2 -|q| 2 )| ( T y | 2 -| ( T S | 2 


Tde 






= |q|((7y<75 + O-^fTy) 


T d f 






= 2|q|g 3 |o-y| 2 


Tdg 


= Tg d 


+ 2|q| 2 |ov| 


2 =-(^M 2 + g 2 -|q| 2 )|^| 2 -|a s | 2 


Tdh 




= T hd 


= -^(°"\>°"5 + o-50- y )M - i(a* v a s - a* s a v )q 3 


Tef 




= Tfe 


= 2^ a v a s ~ (y s a v )M - (a v a s + cr* s a v )q 3 


T 




= T 


= -|q| {^v a s + cr* s a v ) 


T e h 




= T he 


= -M\q\\a v \ 2 


Tfg 




= Tgf 


= -2|q|g 3 |cTy| 2 


Tfh 




= Thf 


= -iHivyVs ~ cr s a v) 


T 

J-ga 






= -Mlqll^yl 2 


Tg h 






= -i|q|(^v-^s - °s°v) 


Tgh 




= T h g 


= 2 ( a v a s + ^*s a v)M + i(aya s - a* s a v )q 3 
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T ha = -\q\(a v a s + a* s a v ) 

Thb = -2|q|g 3 |^| 2 (B.8) 
and for the axivector case the nonzero elements are 

Tab = -T ba = ^((Tya s + <j* s a v )M + i(a v a s - cr* s a v )q 3 

T ac = -T ca = -i|q|((7y(7s - (T* s <7y) 

T ad =-T da = M\q,\\a v \ 2 

T bc =T cb = 2g 3 |q||oy| 2 

T bd =T db =-H{?v°s + o*s°v) 

T c d = -T dc = ^v a s - cr* s a v )M - {a v a s + cr* s a v )q 3 

T ce = -M|q||<7 y | 2 

T c f = -\q\(<Tv<rs + o*s°v) 

T cg = -T gc = --{crycrs ~ a* s a v )M + (a v a s + 0soy)g& 

T ch = T h e - 2|q| 2 Wv\ 2 = -{\m 2 + qi + |q| 2 )M 2 + Ws\ 2 

Tde = ^|q| - 0"5°"v) 

Tdf = 2|q|g 3 |oy| 2 

T dg = T gd + 2|q| 2 \a v \ 2 = -(^M 2 + q 2 - |q| 2 )|a y | 2 + \a s \ 2 

T d h = -T h d = ^v a s - o-* s a v )M - (a v a s + cr* s a v )q 3 

T e f = -T fe = ^(a v a s + a* s a v )M + i(a v a s - cr* s cr v )q 3 

T eg = -T ge = i|q| (a v a s ~ cW) 

T eh =T he = M\q\\a v \ 2 

Tfg =T gf = -2\ci\q 3 \ay\ 2 

T fh = -T hf = -|q|«cr s + a* a v ) 

Tga =-M\ci\\ay\ 2 

Tg b = - 1 q|((Jy a s + a* s a v ) 

1 

Tg h = -T hg = --{<Jy<j s ~ cr* s a v )M + {<Jy(Ts + crsa v )q 3 

T ha = -i\q\(cr v (Js ~ <?* s <?v)i 

Thb = -2|q|g 3 |^| 2 - (B.9) 
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TABLES 



TABLE I. Summary of all possible scalar and vector states in the nonrelativistic limit according 
to the generalised Pauli exclusion principle. 



J 


L 


S 


jAC 


State 











o-+ 


Axiscalar 




1 


1 


0++ 


Scalar 


1 





1 


1— 


Vector 




1 





1+- 


Axivector 






1 


1++ 


Axivector 




2 


1 


1"" 


Vector 
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TABLE II. Bound state masses M (±0.001) for fermion masses from to 5. 



m 


Vector C = +1 

J AC = 1-+ 


Vertor C = — 1 

1 


Axivprtor C = -1-1 
1++ 


Avivprtor C = — 1 
1+- 





0.104 


0.074 


0.119 


0.121 


0.001 


0.112 


0.081 


0.128 


0.129 


0.004 


0.134 


0.099 


0.152 


0.155 


0.009 


0.160 


0.122 


0.179 


0.182 


0.016 


0.194 


0.151 


0.217 


0.223 


0.025 


0.239 


0.184 


0.260 


0.282 


0.036 


0.281 


0.219 


0.295 




0.049 


0.326 


0.258 


0.331 




0.064 


0.372 


0.299 


0.369 




0.081 


0.418 


0.343 


0.410 




0.1 


0.465 


0.389 


0.454 




0.5 


1.341 


1.252 


1.316 




1 


2.325 


2.236 


2.294 




2 


4.344 


4.236 


4.301 




3 


6.357 


6.235 


6.305 




4 


8.385 


8.253 


8.326 




5 


10.370 


10.232 


10.306 
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TABLE III. Comparison between the calculated spectrum of 4-component QED3 positronium 
states in the chiral limit m = and low lying scalar and vector meson masses from Ref. [Oj. The 
h QED3 state is a Goldstone boson generated by the breaking of a global U (2) symmetry. 



jAC Qr jPC 


Positronium mass M 


Meson 


o-+ 





vr°(135) 


0++ 


0.077 


/ (400-1200) 


0~ 


0.111 




0~ 


0.123 




1— 


0.074 


p(770), cu(782) 


1-+ 


0.104 




1++ 


0.119 


ai(1260), /i(1285) 


1+- 


0.121 


/ii(1170), 6i(1235) 



TABLE IV. Charge parities of the Dirac coefficients in the vector, pseudovector, axivector 
and axipseudovector vertex functions. The coefficient functions f n are related to the functions 
a, b, c, d, e, f, g and h by Eq. ( [B.2D and ( |B.3| ). 





C 


Even functions of q.P 


Odd functions of q.P 


Vector 


+1 




hi hi hi hi hi h 




-1 


hi hi hi hi hi h 


hi h 


Pseudovector 


+1 


hi h 


hi hi hi hi hi h 




-1 


hi hi hi hi hi h 


hi h 


Axivector 


+1 


hi hi h, h 


hi hi hi h 




-1 


hi hi hi h 


hi hi hi h 


Axipseudovector 


+1 


hi hi hi h 


hi hi hi h 




-1 


hi hi hi h 


hi hi hi h 
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FIGURE CAPTIONS. 



Figure 1: Bound state masses M against fermion mass m = to 0.1. Full BSE solutions 
from Table|: J AC = 1~+ (O), 1~ (□), l+~ (+) and 1++ (x) states. 



Figure 2: Bound state masses M—2m for fermion mass m = to 1. Full BSE solutions from 
Table|: 1_+ (° solid), 1~ (□ solid), 1+" (+ solid) and 1 ++ (x solid). Nonrelativistic 
predictions of Eq. fl4.2|) using Eq. (|4.3|): (Dashed curves from bottom to top) A' = 
1.7969 (£ = 0), 2.5666 (£ = 1), 2.9316 (I = 0) and 3.1148 (£ = 2) respectively. 
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